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Abstract 

The ARGO-YBJ experiment has been designed to study the Extensive Air Show- 
ers with an energy threshold lower than that of the existing arrays by exploiting 
the high altitude location (4300 m a.s.l. in Tibet, P.R. China) and the full ground 
plane coverage. The lower energy limit of the detector (E ~ 1 GeV) is reached by 
the scaler mode technique, i.e. recording the counting rate at fixed time intervals. 
At these energies, transient signals due to local (e.g. Forbush Decreases) and cosmo- 
logical (e.g. Gamma Ray Bursts) phenomena are expected as a significant variation 
of the counting rate compared to the background. In this paper the performance of 
the ARGO-YBJ detector operating in scaler mode is described and discussed. 

Key words: Cosmic Rays, Gamma Ray Bursts, Resistive Plate Chambers, Air 
Shower Arrays, Extensive Air Showers 
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1 Introduction 



The study of gamma ray sources in tlie 1 - 100 GeV energy range lias so far 
been only partially covered by observations by the satellite-based EGRET, 
with a maximum detectable energy of 30 GeV. On the other hand, present 
ground-based detectors such as Cherenkov Telescopes are still trying to reduce 
their threshold to energies well lower than 100 GeV. 

In this energy range. Gamma Ray Bursts (GRBs) are certainly the most in- 
triguing transient phenomena. GRBs have been deeply studied in the keV- 
MeV energy range, but only EGRET, in the last decade, provided some in- 
formation in the GeV range. Three bursts have been detected at energies 
>1 GeV, with one photon of E=18 GeV [1,2]. The study of the high energy 
emission of GRBs could provide extremely useful constraints to the emission 
models and the parameters of the surrounding medium. 
The interest for transient emission in the GeV region is not only restricted 
to extragalactic phenomena such as GRBs but also includes events closer to 
us - inside our Solar System - such as gamma rays from Solar Flares or high 
energy solar protons producing short time increases of the detector counting 
rate at the ground (Ground Level Enhancements). 

The study of these transient phenomena can be successfully performed by 
ground-based experiments such as Extensive Air Shower (EAS) detectors 
working in "single particle mode" [3] . EAS arrays usually detect air showers 
generated by cosmic rays of energy E > (1 — 100) TeV; the arrival direction 
of the primary particle is measured by the time delay of the shower front in 
different counters, and the primary energy is evaluated by the number of sec- 
ondary particles detected. 
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EAS arrays can work in the energy region E <1 TeV operating in single particle 
mode, i.e. counting all the particles hitting the individual detectors, indepen- 
dently of whether they belong to a large shower or they are the lone survivors 
of small showers. Because of the cosmic ray spectrum steepness, most of the 
events detected with this technique are in fact due to sohtary particles from 
small showers generated by 1 — 100 GeV cosmic rays. 

Working in single particle mode, an air shower array could in principle detect 
a transient emission in the interval 1 — 100 GeV (both by gamma rays or 
charged primaries) if the secondary particles generated by the primaries give 
a statistically significant excess of events over the background due to cosmic 
rays. 

The power of this technique is in its extreme simplicity: it is sufficient to 
count all the particles hitting the detectors during fixed time intervals (i.e. 
every second or more frequently, depending on the desired time resolution) 
and to study the counting rate behaviour versus time, searching for significant 
increases. The observation of an excess in coincidence with a ORB detected 
by a satellite would be an unambiguous signature of the nature of the signal. 
The single particle technique does not allow one to measure energy and di- 
rection of the primary gamma rays, because the number of detected particles 
(often only one per shower) is too small to reconstruct the shower parameters. 
However, it can allow to study the temporal behaviour of the high energy emis- 
sion, and, with some assumptions on the spectral slope (possibly supported by 
satellite measurements at lower energies), can give an evaluation of the total 
power emitted. 

The background rate depends on the altitude and, to a lesser extent, on the 
geomagnetic latitude; it is not strictly constant and several effects both en- 
vironmental and instrumental are responsible for variations on different time 
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scales with amplitudes up to a few per cent. For this technique, an accurate 
knowledge of the detector and its sensitivity to both environmental and in- 
strumental parameters is of crucial importance in order to correctly evaluate 
the statistical significance of the detected signals. 

Obviously the sensitivity of this technique increases with the detection area 
and the energy threshold decreases at higher observation levels, so the most 
suitable detectors have a large area and work at very high altitude. 
The ARGO-YBJ experiment, located at the YangBaJing High Altitude Cos- 
mic Ray Laboratory (Tibet, P.R. China, 4300 m a.s.l.) with a detection area 
of ~6700 m'^, successfully exploits the full coverage approach at very high 
altitude with the aim of studying the cosmic radiation with a low energy 
threshold. 

In order to cover many physics items in different energy ranges, two main op- 
eration modes have been designed: (1) shower mode, with an energy threshold 
of a few hundreds of GeV, and (2) scaler mode, with a threshold of a few GeV, 
this latter differing from the single particle mode for different particle multi- 
plicities (i.e. not only the single particles) are put in coincidence and counted 
on different scaler channels. 

The large field of view (>2 sr, limited only by atmospheric absorption) and 
the high duty cycle (close to 100%) typical of EAS detectors, make ARGO- 
YBJ suitable to search for unknown sources and unpredictable events such 
as GRBs or other transient phenomena. The site location (latitude 30° 06' 
38" N, longitude 90° 31' 50" E) allows sky monitoring in the declination band 
-10° <5 < 70°. 

In this paper the performance of ARGO-YBJ working in scaler mode will be 
described, with some examples of application of this technique in solar physics 
and in the search for GeV emission from GRBs. 



2 The Detector 



The ARGO-YB J experiment is made by a single layer of Resistive Plate Cham- 
bers (RPCs) [4] housed in a large building (100x110 m'^). The detector has 
a modular structure: the basic module is a cluster (5.7x7.6 m^), made of 12 
RPCs (2.850x1.225 m'^ each). 130 of these clusters are organized in a full 
coverage carpet of 5600 m^ with active area ~93%; this central detector is 
surrounded by 23 additional clusters with a coverage of ~40% ( "guard ring" ) 
to improve the core location reconstruction. The installation of the central 
carpet was completed in June 2006. 

Each RPC is read via 80 strips (6.75x61.80 cm^), logically organized in 10 
pads of 55.6x61.8 cm'^ that are individually recorded and that represent the 
space-time pixels of the detector. The RPC carpet is connected to two dif- 
ferent data acquisition systems, working independently, and corresponding to 
the two operation modes, shower and scaler. In shower mode, for each event 
the location and timing of every detected particle is recorded, allowing the 
lateral distribution and arrival direction reconstruction. In scaler mode the 
total counts are measured every 0.5 s: for each cluster, the signal coming from 
the 120 pads is added up and put in coincidence in a narrow time window 
(150 ns), giving the counting rates of >1, >2, >3, >4 pads, that are read by 
four independent scaler channels. The corresponding measured counting rates 
are ~ 40 kHz, ~ 2 kHz, ~ 300 Hz and ~ 120 Hz. 

In order to correctly handle the data, it is very important to evaluate the re- 
sponse to particles hitting the detector. For scaler mode operations the most 
important parameters are the detector efficiency and the strip cross-talk, i.e., 
the probability of having adjacent pads fired by a single particle, giving false 



coincident counts. These effects strongly depend on the adopted gas mixture 
and power supply (in our case Argon 15%, Isobutane 10% and Tetraffuo- 
roethane 75% at 7200V). In these conditions a detector efficiency e = 95% 
has been measured together with a strip occupancy, i.e., the mean number of 
strips fired by a single particle, of 1.2 [4]. A very simple analytical calculation 
has been performed to evaluate the contribution of single particles to higher 
multiplicities: for the given strip occupancy, there is a probability of 10% for 
each strip to generate double counts. A "cross-talk region" , approximated for 
sake of simplicity to a step function, is defined as a width corresponding to 
10% of the strip span. Taking into account the measured detection efficiency, 
the possibility of having cross-talk on both short and long strip sides, and 
the front-end digital logic - giving a maximum number of counts for each pad 
equal to 1 independently of the number of particles hitting simultaneously the 
pad - we obtain the following probabilities of measuring 0, 1 or 2 pads from 
a single particle hitting the RPC: p(0) = 5%; p(l) = 93%; p(2) = 2% (our 
analytical approximation does not give multiplicities greater than 2). This re- 
sult has been checked experimentally on-site with one chamber, finding: p(0) 
= 5.8%; p(l) = 92.0%; p(2) = 2.1%; p(>3) = 0.04% [5]. The good agreement 
confirms the validity of our approach. 

Finally, a substantially different behaviour is expected by the single particle 
channel (> 1) with respect to the higher multiplicity channels, because of the 
contributions given by spurious dark counts (related to front-end electron- 
ics and chamber noise) and natural radioactivity. Since from both we expect 
mostly single counts, these contributions should infiuence only the >1 scaler 
due to the very low rate of random coincidences. An on-site measurement to 
disentangle all these contributions to the single particle rate has been sched- 
uled. 
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3 Detector Stability 



Several effects are known to modify the background counting rate due to 
cosmic rays. First of all, variations in the atmospheric pressure affect the sec- 
ondary particle flux: an increase (decrease) of the ground level atmospheric 
pressure results in a reduction (enhancement) of the background rate, because 
of the larger (smaller) absorption of the electromagnetic component. By mon- 
itoring the pressure at the detector position it is possible to correct the data 
for this effect. Local phenomena connected to the temperature and/or to the 
wind [6] at the level of the early shower development have also been measured, 
the effect being a change of the probability of secondary pions to decay with 
respect to their interaction, due to the local density variation. These effects 
are expected to be negligible compared to the above mentioned barometric 
effect, and very difficult to measure since the outdoor temperature, that is 
usually monitored, is weakly correlated to the high atmosphere (~ 10 km) 
temperature. Moreover, the experiment hall temperature modifies the detec- 
tor performance, changing the gas density inside the RPCs and the bakelite 
electrodes resistivity. 

The cosmic ray background variations due to these parameters are not negli- 
gible (up to a few per cent of the counting rate in some cases), nevertheless 
for signals shorter than a few minutes (such as a high energy emission from a 
GRB) they can be neglected due to the very different time scales of the two 
phenomena. Much more troublesome are possible instrumental effects, such as 
electric noise, that could simulate a GRB signal producing spurious increases 
in the background rate. Working in single particle mode requires very stable 
detectors, and a very careful and continuous monitoring of the experimental 
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conditions. By comparing the counting rates of the scaler channels and re- 
quiring simultaneous and consistent variations in all of these, it is possible 
to identify and reject most of the noise events due to instrumental effects. In 
the following subsection the study of the detector behaviour over short time 
periods is discussed. 

On the other hand, when the study involves long duration signals (as for solar 
activity or sidereal anisotropics) a correction to the counting rates is required 
in order to feature the underlying phenomena. The effects of our corrections 
for the atmospheric pressure and the hall temperature and their efficiency in 
shrinking the counting rate distributions towards a Poissonian distribution are 
discussed in section 3.2. 

3. 1 Detector stability over short time periods 

We define as "short" a time period with a duration less than half an hour. 
Figs. 1 and 2 show the counting rates for multiplicities >1, >2, >3, >4 and 1, 
2, 3 for a typical cluster during 30 minutes of data acquisition. The counting 
rates Cj are obtained from the measured counting rates C>j using the relation: 

a = C7>.-C>,+i (^ = 1,2,3) (1) 

The distributions for all C>j and Cj are Poissonian, meaning that the envi- 
ronmental parameter variations have negligible effects over such a time scale. 
A useful variable is obtained by adding up counting rates of different multi- 
plicities for the same cluster or different clusters for the same counting rate 
multiplicity. In this case, non-Poissonian behaviours appear due to the corre- 
lation between different channels. For the total number of particles hitting a 
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cluster we have: 

4 

Ctot = Y.^C,-Ci + 2C2 + 3C3 + 4C>4 = C>i + C>2 + C>3 + C>4 (2) 

Applying the error propagation to eq. (2) we obtain: 

al^^^ = (Ci) + MC) + 9{C,) + 16(C>4) (3) 

that is larger than a Poisson distribution. Fig. 3 shows the experimental distri- 
bution of the counting rate Ctot for the same cluster and time interval shown in 
figures 1 and 2. As expected, the distribution is Gaussian, wider than Poisso- 
nian and with a width o"exp= (223±4) Hz, compatible with the value of (220±6) 
Hz obtained by applying eq. (3). 

More difficult to estimate is the correlation between different clusters due to 
the probability that some of the counts are generated by the same events giv- 
ing the same multiplicity in different clusters. These events mean the counts 
in different clusters are not independent, widening the distribution from Pois- 
sonian with the sum of counts from more clusters. An analytical calculation 
gives the following result: if Mj is the counting rate of cluster i for a given 
multiphcity (>1, >2, >3, >4 or 1, 2, 3) and n is the total number of clusters, 
we have: 



i=l 

^k* = (^*o*)fl + E(^'-^)p.l (4) 



i=2 



where pi is the probability that the multiplicity is satisfied for the same event 
on (and only on) i clusters. As an example, if we consider three clusters and 
multiplicity >1, we can have two or three clusters hit by at least one particle 
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coming from the same shower. In this case, ps = j^ and p2 = ^'L2+Ni?+N2? ^ 
where iVi23 is the number of counts generated by events hitting all the three 
clusters (with at least one particle) and Nij by events hitting one couple of 
clusters. N^j is not included in A'^123: if three clusters are hit by the same event 
this contributes only to iVi23. 

The width of the Mtot distribution is larger than the Poissonian value crj^^^^ = 
(Mtot)', the Pi values depend on the given multiplicity, the EAS lateral distri- 
bution function and the cluster layout. Since the correlation between different 
clusters is topologically-dependent, for the same number of clusters and mul- 
tiplicity we expect a wider distribution for close clusters compared to more 
distant ones. This has been checked experimentally. As the sensitivity of the 
scaler mode technique is inversely proportional to background fluctuations, 
a continuous carpet is at a disadvantage with respect to scattered detectors 
with the same total surface area. 

The behaviour is the same for the distribution of the ratio between signal and 
background fluctuations (instead of the count rates), as in the GRB search 
analysis. Considering a number s of counts in a time interval At corresponding 
to the GRB duration and a background b of counts averaged on a time interval 
10 times At before and after the signal, we define the normalized fluctuation 
function: 



f = {s-b)/a, (7=^6 + 6/20 (5) 

which represents the significance of the excess compared to background fluc- 
tuations. 

Fig. 4 shows the / distribution for the same cluster of Figs. 1-3 and for At 
= 10 s, obtained using one day of data, for multiplicity > 1 and for the sum 
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of the four multiplicities from > 1 to > 4. Notice that for the single channel 
the Poissonian behaviour of the / distribution is preserved even if the total 
distribution is obtained with 24 hours of data, provided that the total time 
duration of signal plus background is less than 30 minutes. Since the value of 
/ is calculated independently for each signal and background time intervals, 
the long term variations can be neglected and the single channel response is 
always "locally" Poissonian. On the contrary, the / distribution of the sum 
J2t=iC>i follows equation (3), and the width of 1.07 can be easily obtained 
by using the experimental counting rates. The fact that the distribution of 
/ is not affected by long term variations is also very useful to select good 
and bad channels, the latter showing very clearly anomalous non-Poissonian 
fluctuations. 



3.2 Detector stability over long time periods and correlation with environ- 
mental parameters 

As already stated at the beginning of this section, several effects are expected 
to modulate the counting rate measured by our RPC detector, affecting both 
the shower development and instrumental response. These effects must be cor- 
rected before any analysis of long term phenomena, i.e. with duration greater 
than ~ 30 minutes. 

Changes in the atmospheric pressure (P) and the RPC temperature (T) are 
expected to dominate the count rate variations. The correlation with the at- 
mospheric pressure has been widely studied, and a typical coefficient of ~ 
0.7%/mbar is found for a very large range of primary energies [7]. At our 
atmospheric depth (606 g cm^^), the vertical flux of soft (e) and hard (/i) 
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components are almost equivalent [8], reducing the dependency of the pres- 
sure correction coefficient on the atmospheric depth (i.e. zenith angle) due to 
the muonic component. The connection between the counting rates and the 
local temperature is mainly due to the anti-correlation of the bakelite elec- 
trode plates resistivity with the temperature [9,10]. In addition, a variation 
in the detector temperature modifies the gas density, causing a change in the 
effective gap width (the gas volume being constant): the final result is a direct 
correlation between the counting rate and the detector temperature, that we 
assume is the temperature of the gas just out of the RPC 
Fig. 5 shows the counting rates for 3 clusters for multiplicities from > 1 to 
> 4 as a function of local solar time for 72 hours of data accumulation, to- 
gether with the P and T behaviour. Besides the atmospheric pressure and the 
detector temperature, the solar diurnal anisotropy [11] is responsible for the 
counting rate modulation, while the differences in efficiency, chamber dark 
current and electronics noise give the offset displacements. 
In order to correct for these effects we performed a fit with a two-dimensional 
function (as a first approximation linear in P, T) given by the following ex- 
pression: 

CiP,T) ~ Co[l - KP - Po)][l + /3(T - To)] (6) 

where Pq and Tq are the reference values of pressure and temperature for 
counting rate Cq, and /i and j3 are respectively the barometric and thermal 
coefficients. For the counting channel > 1 we found a barometric coefficient 
/i = 0.3-0.5%/mbar and a thermal coefficient f3 = 0.2-0.4%/°C, depending on 
the cluster considered and the specific experimental conditions. For the other 
scalers we obtain ^ = 0.9-1. 2%/mbar and f3 = 0.2-0.4%/°C [12]. 
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Fig. 6 shows the counting rate distribution for C>i and C>2 before and af- 
ter corrections for the previous parameters in a period of 24 hours. It can be 
seen that for the counting channel > 1 the corrected data becomes Gaussian 
but not yet Poissonian, with an experimental standard deviation 23% greater 
than the expected value. In contrast, the distribution of C>2 is Gaussian for 
the raw data and becomes Poissonian after corrections. The fact that the two 
raw distributions for the same day and the /i, j3 coefficients are completely 
different for C>i and C>2 may be explained assuming that part of the C>i 
counts are not due to cosmic rays. 

The statistical behaviour of the counting rates after corrections for atmo- 
spheric pressure and RPC temperature has been checked for longer periods 
showing an agreement better than 5% up to 1 month for C>2 or higher scalers, 
while for C>i the distribution becomes asymmetrical and much larger (~36%) 
than Poissonian already after 48 hours. 

3.2.1 The Forbush Decrease of January 2005 

An example of detection of a long-term phenomenon, linked to solar activity, 
is given by the Forbush Decrease (FD) observed in coincidence with the YBJ 
neutron monitor in January 2005. A FD is the fall in the cosmic ray intensity 
at the Earth one or more days after an Earth directed Coronal Mass Ejection 
(CME) lift-off, and is caused by the shielding effect of the magnetic field and 
shock associated with the high speed plasma cloud emitted from the Sun at 
the time of the eruption. Large solar fiares can be associated with Ground 
Level Enhancements (GLEs) and, depending on the fiare location, may also 
produce a CME that causes a FD; both of these phenomena are normally de- 
tected at the ground by neutron monitors, but a large area detector working 
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in scaler mode is also suitable. 

The FD mentioned above was recorded by the worldwide network of neu- 
tron monitors as a significant decrease in the counting rate after a sudden 
storm commencement. The IGY neutron monitor at Jungfraujoch (Switzer- 
land, vertical effective cut-off rigidity -Rc=4.5 GV) observed a decrease with 
onset around noon UT on January 17, 2005, and with a maximum amplitude of 
about -15% [13]. On January 20, during the recovery phase of the FD at Earth, 
an X7.1 solar fiare created a very intense burst of energetic paticles and a GLE 
was observed, with onset time between 6:48 and 6:57 UT and maximum ampli- 
tude > 5000% at the South Pole neutron monitor station. The GLE spectrum 
was very soft, with a spectral index approaching asymptotically a; ~ -5 [14], 
producing quite small responses at even moderate geomagnetic cut-off rigidity 
locations like Jungfraujoch, where the maximum amplitude was 11.4%. FDs 
similarly produce an energy dependent drop in intensity deeper at low cut- 
off rigidity locations, but since their spectra are much harder than in GLEs 
they should be observable by the ARGO-YBJ array, characterized by Rc=14: 
GV [15]. Fig. 7 shows the YBJ neutron monitor data compared with the four 
counting rate distributions corrected only for the atmospheric pressure since 
the RFC temperature measurement was not yet implemented. Only 12 clus- 
ters were in operation at that time. Even if the scaler mode counting rate is 
more sensitive to environmental parameters, a correlation can be clearly seen, 
with an amplitude of -7% in the YBJ neutron monitor data and about -5% 
and -4% for C>i and C>2 data. No decrease is observed for C>3 and C>4, and 
no GLE is seen in any scaler channel and in the YBJ neutron monitor data. 
Another interesting feature of FDs is the recovery phase. While the magnitude 
of FDs depends on the cut-off rigidity (as already shown), the recovery time is 
expected to be independent of the galactic cosmic ray energy and thus similar 
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for all detectors. Despite of this, as pointed out in [16], the recovery time in 
some events strongly depends on the median energy of the detected particles, 
with a faster recovery for more energetic cosmic rays. Moreover, these higher 
energy cosmic rays show an " over- recovery" , i.e. a recovery of the cosmic ray 
intensity to a level higher than that before the FD. The study of these features 
is very important to characterize and model the propagation of the interplan- 
etary shocks; due to its high cut-off rigidity, ARGO-YBJ can contribute to 
the study of acceleration processes, keeping in mind that the energy upper 
limit for incoming solar particles in the Earth environment has not yet been 
established. Due to its huge area, this study can be done with high sensitivity 
and small integration time. For example, a distinct enhancement of the diur- 
nal variation in the recovery phase after the FD can be seen mainly around 
24-25 January. It is most obvious in the C>2 data, but it is also evident in the 
C>3 and C>4 data where the FD is not observed. This phenomenon is known 
but not extensively studied, hence our data can give valuable information on 
the energy dependency of the FD recovery. A more detailed analysis of these 
data and their correlation is in progress. 



4 Effective Area and Expected Counting Rates 

The effective area is usually obtained by simulating a large number of showers 
^sim with the core sampled on a wide area Agim-, much wider than typical 
shower dimensions, and counting the number of showers Ntrig matching the 
trigger condition: 

^eff{E) = Asim ■ Ntrig /N Sim (7) 



For the ARGO-YBJ detector, the effective area has been estimated by means 
of a detailed Monte Carlo (MC) simulation for photons, protons and helium 
of fixed energies ranging from 10 MeV to 100 TeV and fixed zenith angle 9 
from to 50 degrees (with 10 degrees steps). The CORSlKA/QGSJet code 
6.204 [17] has been used with a full electromagnetic component development 
down to Ej/jr= 1 MeV for electrons and photons and 50 MeV for muons and 
hadrons. 

Since the actual efficiency depends essentially on the shower particle lateral 
distribution, a huge number of showers must be simulated over an area large 
enough to completely contain the event. 

In order to save computing time the shower sampling has been performed by 
means of the "reciprocity technique" [18]: the sampling area A5 is uniformly 
filled with replicas of the same ARGO-YBJ carpet, one next to the other, 
for a total area of ~ 10"^ x 10^ m^. Following the reciprocity technique, we 
sample the shower axis only over the area covered by the carpet located at 
the center of the array, with the prescription of considering the response of 
all the detector replicas. Then, event by event, we calculate the number of 
clusters on the entire grid which contain at least 1, 2, 3, 4 fired pads. In this 
calculation the cross-talk between adjacent pads has been taken into account, 
as discussed in section 2. The resulting effective areas for the full ARGO-YBJ 
detector as a function of the energy of primary photons and protons at zenith 
angle 9 = 20° are shown in Fig. 8 in the four multiplicity channels (1, 2, 3 and 
>4). 

If the probability that more than one particle per shower hits the detector 
is negligible, i.e. for sufficiently small detector and/or small primary energy, 
the effective area Agcaier scales like the mean electron size after the shower 
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maximum at the detection level [3]: 

Ascaler ~ 6 ■ {Nch) ' Ad ■ COS 6. (8) 

where e and A^ are the detector efficiency and area, and {Nch) is the mean 
number of charged secondary particles. /,From Fig. 8 we see that, as expected, 
the dependence is described by a power law index ^1.6 for multiplicity 1, and 
this index increases with the multiplicity approaching the value ~ 2.6, valid 
only in the shower trigger configuration [19]. The index lowers with increasing 
energy when the formula (8) loses its validity because the area of the cluster is 
too large for the given shower particle density and thus more than one particle 
ffies it. The effective area for protons is usually lower than that for photons of 
the same energy, increasing the sensitivity of the measurement, except in the 
10 - 100 GeV energy range, where, due essentially to the muonic component, 
the total number of charged particles in proton-induced showers is larger than 
that in photon events. Moreover, the geomagnetic cutoff at YBJ is >14 GV 
(depending on the arrival direction) [15]. 

In order to get the counting rates, the effective area is folded with the primary 
spectrum, expressed by a differential power law with spectral index a = -2.0 
for photons, a = -2.77 for protons and a = -2.64 for helium [20]. We adopted 
a = -2.0, which is the mean spectral index for GRBs as measured by EGRET 
[21]. The resulting convolutions for photons and protons are shown in Fig. 9. 
The behaviour of these curves shows how each energy decade contributes to 
the counting rate for a fixed zenith angle. One can note that using four scalers 
instead of a single one allows the investigation of different energy ranges, with 
the possibility of giving an indication on the source spectral index in case of 
positive detection. 
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The resulting counting rates, integrating the different effective areas for an- 
gular bins up to a zenith angle of 60 degrees and both proton and helium 
contributions, are 24.4 kHz for Ci, 1.8 kHz for C2, 195 Hz for C3 and 117 Hz 
for C>4. The uncertainties on these rates come mainly from those affecting the 
primary fluxes which are estimated to be about 20% (see for example [22]). 
A comparison with the measured rates, i.e., 38 kHz for Ci, 1.7 kHz for C2, 
180 Hz for C3 and 120 Hz for C>4, shows a very good agreement between 
the predicted values and the experimental, apart from the counting rate Ci 
affected by spurious signals as explained in section 2. 

5 Upper Limits and Expected Sensitivity 

The energy fiuence F in the energy range [-Emm, E^ax] corresponding to an 
excess of N events above the cosmic ray background is given by: 

F(a, Emax) = -p — -^^ ; ^ 9 

SE:::E-^eff{E,e)dE 

where Aeff{E, 6) is the effective area calculated in a given multiplicity channel 

at the fixed zenith angle 6 and a is the photon spectral index (see previous 

section). 

In case of no significant excess, typical fiuence upper limits to GRB emissions 

are obtained using the measured counting rates in each of the four multiplicity 

channels and a model for the GRBs. 

Assuming a power law spectrum up to Emax ranging from 100 MeV to 1 

TeV, no 77 absorption by the Extragalactic Background Light, a zenith angle 

6 = 20° and a duration At = 10 s, the upper limits corresponding to a Acr 
excess in the multiplicity channel 1 are calculated as a function of the spectral 
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index in the energy range 10 MeV — Emax- The results are reported in Figure 



10. For different time durations, upper hmits scale as w At(s)/10. 
In order to evaluate the real possibilities of a GRB detection, we compare the 
ARGO-YBJ sensitivity with the GRB fluxes measured by EGRET. During 
its lifetime EGRET detected 16 events at energies >200 MeV[l]. Above this 
energy, all their spectra show a power law behaviour without any cutoff, sug- 
gesting the possibility that a large part of GRBs could emit GeV or even TeV 
gamma rays. 

In our simple model every EGRET GRB spectrum is extrapolated to higher 
energies, using the measured power law index, up to a maximum energy Emax 
representing a possible intrinsic cutoff at the source. Then we consider an ex- 
ponential cutoff to take into account the effects of the extragalactic absorption, 
obtaining the spectrum: dN/dE = K E" e~'^^^'^\ The intrinsic energy cutoff 
is assumed in the range 100 GeV < E^ax < 1 TeV and the GRB redshift in 
the range 0< z <3. The absorption factor is calculated by using the values of 
optical depth t{E,z) given in [23]. 

Our calculation shows that a fraction ranging from 10% to 50% of EGRET 
GRBs is above the ARGO-YBJ minimum detectable flux, depending on the 
cutoff energy Emax and the redshift z of the GRB. Assuming Emax=^ TeV 
(-Emaa;=100 GcV) and no absorption of gamma rays, ARGO-YBJ could detect 
50% (30%) of EGRET GRBs. If z=l (z=3) this value reduces to ~ 20% (10%), 
for both values of Emax- Considering that EGRET observed 16 GRBs during 
7 years (1991-1997) and that the EGRET fleld of view (~ 0.5 sr) is about 
four times smaller than that of the present Swift satellite, ~ 9 GRBs similar 
to those seen by EGRET are expected to occur in one year in the Swift fleld 
of view. 
Since about 1/10 of these bursts are, by chance, in the ARGO-YBJ fleld of 
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view, we expect ~ 1 EGRET-like GRB detected by Swift to occur in our de- 
tector field of view in one year. A fraction of them (depending on z and Emax) 
will be effectively detectable. 



6 Conclusions 

The ARGO-YBJ detector operated in scaler mode has a high duty cycle and 
optimal stability, fulfilling the fundamental requirements of the scaler mode 
technique. 

The main goal is the measurement of the high energy tail of GRB spectra, 
and the detection of few events during the experiment life time is expected 
if the spectrum of the most energetic GRBs extends at least up to 100 GeV. 
In this case, with four measurement channels sensitive to different energies, 
valuable information on the high energy spectrum slope and possible cutoff 
may be obtained. 

In any case, even if no signal is detected, the large field of view and the high 
duty cycle, allowing for continuous monitoring overhead, make the ARGO- 
YBJ experiment one of the most sensitive detectors for the study of the 
high energy spectrum of GRBs, with typical fluence upper limits down to 
lO^^erg/cm^ in the 1 - 100 GeV region, well below the energy range explored 
by the present generation of Cherenkov telescopes. 

On the other hand, due to the high rigidity cut-off, the large area of the de- 
tector and its high mountain location allow the study of the solar modulation 
of galactic cosmic rays and energetic relativistic particles coming from the 
Sun with high sensitivity and fine time resolution. This has been shown by 
the observation of the January 2005 FD with a preliminary discussion on the 
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recovery phase. 

Data collected up to May 2008 have been analyzed according to the proce- 
dures outlined in this work and the results of these analyses will be presented 
in later papers. 
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8 Captions list 



Fig. 1 Experimental distributions of the counting rates and their Gaussian fits 
for a typical cluster: (a) C>i, (b) C>2, (c) C>3 and (d) C>4 for 30 minutes data 
accumulation. The standard deviation of the Gaussian fit {(Texp) is compared 
with the square root of the mean of the experimental distribution {(Jth) to 
check the compatibility with the Poisson distribution. 

Fig. 2 Experimental distributions of the counting rates and their Gaussian fits 
for a typical cluster: (a) Ci, (b) C2 and (c) C3 for 30 minutes data accumula- 
tion. The standard deviation of the Gaussian fit {(Xexp) is compared with the 
square root of the mean of the experimental distribution (ct/i) to check the 
compatibility with the Poisson distribution. 

Fig. 3 Experimental distribution of the total number of particles hitting a 
typical cluster and its Gaussian fit for 30 minutes data accumulation. The 
standard deviation of the Gaussian fit {(Jexp) is compared with the expected 
value derived from the single multiplicities counting rates (see text, eq. (3)). 
Fig. 4 Experimental distribution of the normalized excesses of signal over back- 
ground of a typical cluster (see text, eq. (5)) compared with a Gaussian fit. 
Top: C>i channel; bottom: sum of the four channels. 

Fig. 5 Counting rates as a function of time for 3 clusters for multiplicities from 
> 1 to > 4 and 72 hours of data accumulation, together with pressure (P) and 
RPC temperature (T) behaviour for the same time interval (lower plots). 
Fig. 6 Experimental distributions of counting rates for C>i and C>2 of a typ- 
ical cluster before (left) and after (right) corrections for pressure and RPC 
temperature in a period of 24 hours. The standard deviation of the Gaussian 
fit {cexp) is compared with the square root of the mean of the experimental 
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distribution (crth) and the resulting "non-Poissonian" behaviour — ^ = '^^"'^ ""*'' 
is reported. 

Fig. 7 Plot of the YBJ neutron monitor (NM) data compared with the four 
counting rate distributions (added on all the operating 12 clusters and cor- 
rected only for the atmospheric pressure) during the Forbush Decrease of 
January 2005 (Neutron Monitor data: courtesy of the Yangbajing Neutron 
Monitor Collaboration). 

Fig. 8 Effective areas for photon and proton-induced showers (left and right 
plot, respectively) with zenith angle 6 = 20° in the four multiplicity channels 
(from top to bottom, 1, 2, 3 and > 4). 

Fig. 9 Convolutions of the effective areas of Fig. 8 with a differential power law 
spectrum of index a = -2.0 for photons and a = -2.77 for protons, taking into 
account in this latter case the geomagnetic bending, in the four multiplicity 
channels (from top to bottom, 1, 2, 3 and > 4). 

Fig. 10 Fluence upper limits in the energy range 10 MeV - E^ax corresponding 
to a 4o" excess in the multiplicity channel 1, as a function of the GRB spectral 
index, for a zenith angle ^=20°, a time duration At = 10 s and a redshift z=0. 
Each set of points refers to different E^ax as shown in the figure. 
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Fig. 1. Experimental distributions of the counting rates and their Gaussian fits 
for a typical cluster: (a) C>i, (b) C>2, (c) C>3 and (d) C>4 for 30 minutes data 
accumulation. The standard deviation of the Gaussian fit (cTf^xp) is compared with 
the square root of the mean of the experimental distribution {(7^^) to check the 
compatibility with the Poisson distribution. 
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Fig. 2. Experimental distributions of the counting rates and their Gaussian fits for 
a typical cluster: (a) Ci, (b) C2 and (c) C3 for 30 minutes data accumulation. The 
standard deviation of the Gaussian fit (cJexp) is compared with the square root of 
the mean of the experimental distribution {(Jth) to check the compatibility with the 
Poisson distribution. 
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Fig. 3. Experimental distribution of the total number of particles hitting a typical 
cluster and its Gaussian fit for 30 minutes data accumulation. The standard devi- 
ation of the Gaussian fit (aexp) is compared with the expected value derived from 
the single multiplicities counting rates (see text, eq. (3)). 
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Fig. 4. Experimental distribution of the normalized excesses of signal over back- 
ground of a typical cluster (see text, eq. (5)) compared with a Gaussian fit. Top: 
C>i channel; bottom: sum of the four channels. 



32 



45000 




18 36 54 

YBJ Solar Time (h) 



18 36 54 

YBJ Solar Time (h) 



72 






180 



200 



18 36 54 

YBJ Solar Time (h) 



o"' 140- 



DC 



o 




18 36 54 

YBJ Solar Time (h) 



Si 

E 




18 36 54 

YBJ Solar Time (h) 



72 



3 



0) 

a. 

E 

0) 

O 
Q. 
DC 




18 36 54 72 

YBJ Solar Time (h) 



Fig. 5. Counting rates as a function of time for 3 clusters for multiplicities from > 1 
to > 4 and 72 hours of data accumulation, together with pressure (P) and RPC 
temperature (T) behaviour for the same time interval (lower plots). 
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Fig. 6. Experimental distributions of counting rates for C>i and C>2 of a typical 
cluster before (left) and after (right) corrections for pressure and RFC temperature 
in a period of 24 hours. The standard deviation of the Gaussian fit {cfexp) is compared 
with the square root of the mean of the experimental distribution {(Jth) and the 
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is reported. 



34 



N 

X 
0) 

ra 
DC 




UT (day of January, 2005) 



X 


490 
480 


O^' 


470 


^sll 


460 




18 20 22 

UT (day of January, 2005) 



^ 24000 f^ 
E, 23500 
J' 23000^ 
[^ 1 225001- 
" 22000 '- 




18 20 22 24 

UT (day of January, 2005) 



I 4200 P 

CO 

Al 

O ^ 4100 
" 4000 




UT (day of January, 2005) 



N 

X 



o 




£ 1740?- 



18 20 22 

UT (day of January, 2005) 



Fig. 7. Plot of the YB J neutron monitor (NM) data compared with the four counting 
rate distributions (added on all the operating 12 clusters and corrected only for 
the atmospheric pressure) during the Forbush Decrease of January 2005 (Neutron 
Monitor data: courtesy of the Yangbajing Neutron Monitor Collaboration). 
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Fig. 8. Effective areas for photon and proton-induced showers (left and right plot, 
respectively) with zenith angle Q = 20° in the four multiplicity channels (from top 
to bottom, 1, 2, 3 and > 4). 
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Fig. 9. Convolutions of the effective areas of Fig. 8 with a differential power law 
spectrum of index a = -2.0 for photons and a = -2.77 for protons, taking into 
account in this latter case the geomagnetic bending, in the four multiplicity channels 
(from top to bottom, 1, 2, 3 and > 4). 
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Fig. 10. Fluence upper limits in the energy range 10 MeV - Emax corresponding to 
a 4(7 excess in the multiphcity channel 1, as a function of the GRB spectral index, 
for a zenith angle 0=20°, a time duration At = 10 s and a redshift z=0. Each set 
of points refers to different E^ax as shown in the figure. 
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